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Abstract 

We calculate the new dinamic exponent 9 of the 4-state Potts model, using short-time simula- 
tions. Our estimates 0\ = —0.0471(33) and 62 = —0.0429(11) obtained by following the behavior of 
the magnetization or measuring the evolution of the time correlation function of the magnetization 
corroborate the conjecture by Okano et. al. In addition, these values agree with previous estimate 
of the same dynamic exponent for the two-dimensional Ising model with three-spin interactions 
in one direction, that is known to belong to the same universality class as the 4-state Potts model. 
The anomalous dimension of initial magnetization xq = z9 + [5/v is calculated by an alternative 
way that mixes two different initial conditions. We have also estimated the values of the static 
exponents (3 and v. They are in complete agreement with the pertinent results of the literature. 
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I. INTRODUCTION 



Several results about critical phenomena have been recently obtained using Monte Carlo 
simulations in short-time regime l|. Such simulations are more convenient than traditional 
ones done in the equilibrium regime because they circumvent a known problem in compu- 
tational physics: the critical slowing down phenomena. 

Investigating Monte Carlo simulations before attaining equilibrium permits to us obtain- 
ing the static critical exponents {(3 and v) but also leads to the less known dynamical ones. 
The scaling equation for non-equilibrium regime was obtained by Jansen et al. [2j], on the 
basis of renormalization group theory. For systems without conserved quantities like energy 
and magnetization (model A in the terminology of Halperin and Hohenberg jjH), it is written 
as 

M (k \t,r,L,m ) = b^M( k \b- z t,b»T,b- l L,b Xo m ) (1) 

r-i 

where mo is the initial magnetization, xo = Xo(mo) m * s the anomalous dimension of the 
initial magnetization, j3 and v are the known static exponents and z is the dynamic one (r ~ 
£*). The values of are the kth moments of magnetization, defined by M' fe ' = (M fc ), 
where (•) indicates an average over several samples randomly initialized but satisfying the 
condition of having the same magnetization (mo) at the beggining. The relevance of this 
result is allowing to include the dependence on the initial conditions of the dynamic systems 
in their non-equilibriun relaxation. As an important consequence, they could advance the 
existence of a new critical exponent 8, which is independent of the known set of static 
exponents and even of the dynamic exponent z. 

By considering large systems and choosing r = (T = T c ) ; b~ z t = 1 in the equation (JTJ), 
we obtain the dynamic scaling law for the magnetization: 

M{t, m ) = t^M(l, t^mo). (2) 

Expanding M around the zero value of the parameter u = t * mo, one obtains: 

M(t,m ) = m t e + 0{u 2 ) (3) 

which leads to the power law: 

M{t) ~ m t e (4) 

_ z 

whereas t < to ~ m x ° since in this regime « C 1. This new universal stage, characterized 
by the exponent 6 = (xo — ) /z , has been exhaustively investigated to confirm theoretical 
predictions and to enlarge our knowledge on phase transitions and critical phenomena. 



2 



Indeed the anomalous behavior of the magnetization at the beginning of relaxation, also 
called critical initial slip, was originally associated to a positive value of 9. This kind of 
behavior was confirmed in the kinetic q = 2 and 3— state Potts models [lj], as well as in 
irreversible models like the majority voter one |5[ and the probabilistic cellular automaton 
proposed by Tome and Drugowich de Felfcio to describe part of the immunological system (f| . 
However, the exponent 9 can also be negative. This possibility was observed in the Blume- 
Capel model (analytical 7] and numerically [^]) and in the Baxter- Wu model (numerically 
P]), an exactly solvable model which shares with the 4— state Potts model the same set 
of critical exponents. In addition, a negative but close to zero value for 9 was obtained 
for the two-dimensional (2-D) Ising model with three-spin interactions in just one direction 
(IMTSI) jljj], a model which is also known to behave to the 4— state Potts model universality 
class. 

However, at the best of our knowledge a numerical determination of the exponent 9 for 
this special case (q = 4) of the Potts model continues to be lacking. Thus, we decided to 
investigate the short-time behavior of that model in order to better understand its dynamics 
and also to learn about the role of the exponent xq of the nonequilibrium magnetization 
concerning the well established classification of the models (at least in what concerns the 
static behavior) in the universality classes. 

As we mentioned above the estimates for the exponent 9 for two models that belong 
to the same universality class as the 4-state Potts model have revealed discrepant results. 
Whereas the two-dimensional Ising model with threee spin interactions in one direction 
IMTSI) obeys the conjecture by Okano et al. |ll|(^ should be negative and close to zero) 



lol | the Baxter- Wu model [9] strongly disagrees from that prediction. Putting in numerical 
values, whereas the IMTSI model exhibits a small and negative value (—0.03 ± 0.01) for 9, 
when studying the BW model we found 9 = —0.186 ± 0.002. 

In this paper we study the short-time behavior of the two-dimensional four-state Potts 
model and present numerical estimates for 9 using two different methods. First, we fix the 
initial value of the magnetization and follow its time evolution to obtain 9(m ) which in turn 
can be extrapolated to lead to 9 when mo — > 0. Second, we deal with the time correlation 
of the magnetization defined by 0|: 

C(t) = ± (M(t)M(O)) (5) 
where the average is done over samples whose initial magnetizations are randomly choosen 
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but obey (M(0)) = 0. In reference |5[ it was shown that this quantity exhibits at T = T c 
the power-law behavior: 

C(t) ~ t 9 . (6) 

This approach has several advantages when compared to the other technique but its appli- 
cation was initially restricted to the models which exhibit up-down symmetry. Recently, 
Tome has shown that this result is more general and can be applied whenever the model 
presents any group of symmetry operations related to the Markovian dynamics. Equation 
|B] was shown to be valid for instance in the antiferromagnet ordering model, models with 



one absent state and eve, in the B aX te,Wn m „de.l which exhibits a Z(2) 3 Z<2) 
symmetry. 

It is well known that the q ^ 2 Potts Model does not have up-down symmetry. However, 
based on the paper of Tome 



12j we can use this approach to calculate the exponent 9 of the 



three- and four-state Potts models. 

We checked this possibility working with the three-state Potts model. Our estimate is 
9 = 0.072(1), in good agreement with the result obtained by Zheng JlJ] (9 = 0.070(2)) and 
with the estimate obtained by Brunstein and Tome 1 1 3] using a cellular automaton which 
exhibits C3v symmetry, the same as the three-state Potts model. A plot of the correlation 
function of the magnetization C(t) in this case is presented in the Fig. 1. 




FIG. 1: Log-log plot of C(t) x t to 3-states Potts Model 

In the sequence we estimated the exponent 9 of the 4-state Potts model which agrees 
with the conjectured by Okano et al. In addition, they agree with previous estimates for 9 
obtained by Simbes and Drugowich de Felicio jicl ]. 



4 



To confirm our first estimate we simulated the magnetization evolution for four different 
initial values of m which after extrapolated (m — > 0) led to a similar result. 

It is clear that having the exponent 9 the anomalous dimension xq follows. But inspired 
in previous studies [14J we decided to investigate a direct manner in determining xq using 
mixed initial conditions. In order to build the adequate function we remember that when 
in contact with a heat bath at T = T c the magnetization of completely ordered samples 
(mo = 1) decays like a power law 

M(t)~t^. (7) 
Thus, according to the scaling relations (JH) and (J7J) it is enough to work with the ratio 



to obtain a power law which decays as t e+l3 / vz = t x °l z . Using z as input [14] we can achieve 
the value of Xq. On the other hand, by following the relation (J7|) we estimate the ratio 
(3/vz which can be compared with the exact result (3/v = 1/12 after using the value of the 
dynamical exponent z. 

Finally using derivatives of the magnetization at an early time and once more the value 
of z as input we obtain the critical exponent v of the correlation length, whose numerical 
estimates by usual techniques (phenomenological renormalization group, hamiltonian studies 
and equlibrium Monte Carlo simulations) are always very different from the pertinent result 

= 2/3). 

The paper is organized as follows: in the next section we present a brief review and some 
details of the simulation. The results are presented in Section 3 and our conclusions are in 
section 4. 



II. THE KINETIC 4-STATES POTTS MODEL 



In time-dependent simulations we are interested in finding power laws for physical quan- 
tities even when the system is far from equilibrium. In this regimen the magnetization M (t) 
must be calculated as an average over several samples because the system does not obey 
any a priori probability distribution. The average can be done in different ways: we can 
prepare all the samples with the same initial magnetization (sharp preparation) or generate 
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samples which satisfy a less restrict criteria like to have mean value of the magnetization 
equal to zero. 

The q states Potts model ferromagnetic without the presence of an external field is defined 
by the Hamiltonian Q: 

H = -jJ2^ (9) 

where J denotes the interaction between the nearest neighbors and a\ can assume 

different values = 0, 1.., q — 1. If two spins are parallel they contribute with energy —J, 
else the energy is null. 

The critical temperature of this model, is known exactly [15] , 

^ = log[l + v^- (10) 

The magnetization, different of other models as Ising and Blume Capel is not only the 
sum of variables of spin. A general expression used for the magnetization in the Potts model 
that considers an average over the sites and over the samples is written as 

< J W» = jy.(g-i)iv gD«W>-'- 1 > (") 

where denotes the spin i of the jth. sample at the tth. MC sweep. Here N s denotes 
the number of the samples and L d is the volume of the system. This kind of simulation is 
performed Nb times to obtain our final estimates as a function of t. In order to prepare a 
lattice with a given magnetization we need to choose the states Oi = 0, 1,2,3 in each site 
with equal probability (1/4). Next we measure the magnetization and change states in sites 
randomly chosen in order to obtain a null value for the magnetization. Finally we change S 
sites occupied by = 0, 2 or 3 and substitute by Oi = 1. The initial magnetization in this 
case will be given by: 

mo = ^ (12) 

where N = L d . 

We have chosen to update the spins according to the heat bath algorithm, which means 
that the probability that the spin a^t localized at the site ~x can assume another value aL 
is given by: 

exp(-^E^ w , 4 ) 

p{*? - 4) = (13) 



E exp(-^E^ Wj4 ) 

aL,=0 k 
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where ^ denotes the sum about the nearest neighbors sites of the spin at the site x . 

~k 

III. RESULTS 

We performed Monte Carlo simulations for Nb = 5 different bins and four different initial 
magnetizations, m = 45 /3N w 0.033, 0.049, 0.066 and 0.082, for a lattice of size L = 90. 
These magnetizations correspond respectively to 5 = 200, 5 = 300, 5 = 400, 5 = 500 
according to the expression (jl2j) . To estimate the errors we ran 5 different bins, with 35000 
samples each one. 




FIG. 2: Time evolution of the magnetization for tuq = 0.033. 

In the figs. 2 and 3 we present the log-log plot for two different initial magnetizations 
and two different values of mo = 0.033 and = 0.082 respectively. The fit had excellent 
goodness (Q = 0.99) in the entire range [U,tf] C [0, 200]. We have chosen the range [10, 100] 
to estimate the value to 9 for all initial magnetizations. We must note that the slope is 
positive for m = 0.082 but changes the signal at some value between 0.066 and 0.049. 
This trend to a negative value can be also certified in the Fig. 4 that shows the plot of 
9 = 9(m ). It is worth to mention that there is a clear linearity in the plot of 9 versus the 
initial magnetization (see fig. 4) but big fluctuations are observed in the evolution of the 
magnetization in each case. The value found for 9 ext is —0.0471(33) is in fair agreement with 
the result Q] for the IMTSI model and strongly disagrees with the result obtained for the 
Baxter- Wu model. 

We also ran Monte Carlo simulations to obtain the time-evolution correlation C(t), and 
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FIG. 3: Time evolution of the magnetization for m,Q = 0.082. 

the the dynamic exponent 9 according to the relation For the same interval [10, 100] we 
have gotten the estimate 9 = —0.0429(11) that corroborates the above estimate —0.0471(33) 
obtained from the evolution of the non-equilibrium magnetization. We used 35000 samples 
and once more 5 different bins to make this estimate. In figure 5 it is shown the log-log plot 
of C(t) x t for the four-state Potts model. 
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FIG. 4: as function of uiq in the 4-state Potts Model (extrapolation) 

Equilibrium simulations for the four-state Potts model are ever acompannied by strong 
fluctuations. The reason is the presence in the Hamiltonian of a marginal operator charac- 
terized by an anomalous dimension equal to zero. If we hope some signal of that anomalous 
behavior in short-time simulations done far from equilibrium we should look for that in the 
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only new anomalous dimension: that of the initial magnetization xq. In order to confirm the 
null value for the anomalous dimension we estimate directly Xq using the reported function 
F 3 . Over here we used = 5 bins for each initial condition and again the N s = 35000 



4-states Potts Model 




t 



FIG. 5: Log-log plot of C(t) X t to 4-states Potts Model 

realizations. A plot when the initial magnetization is m = 0.082 can be seen in the fig. 
6. Extrapolating to m — > we found {xo/z) ext = 0.0077(33) with the same set initial 
magnetizations mo that was used to extrapolating m . Using as input z = 2.290(3), we 
obtain xo = 0.0176(75). 

The value of (3/v can be obtained directly of magnetization decay of a initial ordered 
state, considering the relation: 

\vzj c 

The value of the exponent obtained from the decay (M ~ t~~^), was 

(JL) c = 0.0547(1) with 

Q = 0.84. Using as the input the same value for z value we achieved ^ = 0.12526(28) that 
is in a good agreement with the exact result ^ = 0.125 Q|. 

Finally to find the exponent u, we calculated the expected decay of derivative 



D(t) = A In M(t,r) 



r=0 

(15) 



In M(t, t + 6) - In M(t,T- 6) 
= 26 ' 

According to the reference (JT|, for mo = 1, we expect of D(t) the power law behavior 

D{t)~t™ (16) 
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FIG. 6: Plot of F3 function for case tuq = 0.066. 



1 \ 1 

and so v = ( — I -.In the Fig. 7 we illustrate the power law (1161). performing simulations 
1 vz } „ z 



with N s = 8000 samples with iV b = 5 to estimate errors and t r 



1000 MC steps. Our 




FIG. 7: Time dependence of derivative D(t) of the 4-state Potts Model. The slope gives 1/uz. 

best estimative is — = 0.65453(46), in the interval [10,1000] with Q = 0.99. The value 

found to v at this interval is v — 0.667(1). The value ( — J also is used to estimate 8. 
This minimizes the errors because: 



- 1 ■) (' 



VZ I \ vz , 
/ c \ / c 



(17) 



The value found to 8 is 0.0836(2). These values are in complete agreement with conjec- 
tured results v — 2/3 = 0.6 and 8 = 1/12 = 0.083. 
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IV. CONCLUSIONS 



We calculated the dynamic exponent 9 of the four-state Potts Model using two different 
techniques: the time-evolution of the magnetization when the samples are prepared with a 
small and nonzero value of mo for four different values of initial magnetization and through of 
the time correlation of the magnetization. The values found for 9 in both cases corroborate 
the conjecture by Okano et al. and are in complete agreement with the estimate obtained 
for another model which is in the same universality class as the 4-state Potts model: the 
2-D Ising model with three-spin interactions in one direction. We also estimated directly the 
anomalous dimension of the initial magnetization x$ which results very closed to zero. The 
same does not occur in the Baxter- Wu model that, although exhibiting the same leading 
exponents as the 4-state Potts model, does not own a marginal operator which prevents us 
from determine with good precision the critical exponents. 

We also estimated the static exponents (3 and v using the short-time Monte Carlo simu- 
lations and our results are in surprisingly agreement with pertinent results. 

As an additional contribution we present a new estimate for the exponent 9 of the 3-state 
Potts Model working with the time evolution of the time correlation of the magnetization 
C (t) . This kind of approach has been previoulsy used by Brunstein and Tome when studying 
a four-state cellular automaton with C-3^ symmetry |l3j].The result agrees with estimates 
obtained from the evolution of the magnetization 1[ after the extrapolation (mo — > 0). 
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